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ABSTRACT 

We present a theoretical investigation of the hght, radius and radial velocity 
variations of the prototype 6 Cephei. We find that the best fit model accounts for 
luminosity and velocity amplitudes with an accuracy better than 0.8a, and for 
the radius amplitude with an accuracy of 1.7a. The chemical composition of this 
model suggests a decrease in both helium (0.26 vs 0.28) and metal (0.01 vs 0.02) 
content in the solar neighborhood. Moreover, distance determinations based on 
the fit of light curves agree at the O.Scr level with the trigonometric parallax 
measured by the Hubble Space Telescope (HST). On the other hand, distance 
determinations based on angular diameter variations, that are independent of 
interstellar extinction and of the p-factor value, indicate an increase of the order 
of 5% in the HST parallax. 

Subject headings: Cepheids - stars: distances - stars: evolution - stars: oscilla- 
tions 
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1. Introduction 



X 



Classical Cepheids are the most important primary distance indicators. One of the 
key issues concerning the use of Cepheids as distance indicators is their dependence on 
chemical composition. Empirical tests seem to suggest that metal-rich Ceph eids are, at fixed 
period, brighter than metal-poor ones, either over the entire period range (IKennicutt et al. 



19981 . l2003l : iKanbur et al.ll2003l : IStorm et al.ll2004J : ICroenewegen et al.ll2004l : ISakai et al. 



2004 
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Macri et al.l 120061) or for periods shorter than ^ 25 days (iSasselov et aLlll997l : iSandage et ah 
2004 ). On the other hand, distance estimates based on the Near- Infrared Surface Brightness 
method indicate that the slope of the P eriod-Luminosity (P L) relation might be the same 



for Magellanic and Galactic Cepheids (JGieren et al.l l2005l). However, the impact of the 



proje ction factor p on the Cepheid distances is still debated (ISaha et al.ll2006l : iNardetto et al. 



1999 



2007 ) . Early theoretical predictions based on lin ear pulsation models (e.g. 



Alibert et al. 1999: Baraffe fc Ahbert 



nonlinear convective models (IBono et al. 



2001 



1999 



Sandage et al. 



predicted a mild metallicity effect, while 



Caputo et all l2002l : iFiorentino et all 12002 



Marconi et al.ll2005l ) showed that the metallicity affects both the zero-point and the slope of 
the PL relations. This effect is more evident in the optical bands and it is nonlinear when 
moving from metal-poor [Z = 0.004) to super metal-rich [Z = 0.04) structures. Current 
predictions also suggest a turnover across solar metallicity (Z ^0. 02) due to the correlated 
increase in He content (IFiorentino et al.l 120021 : iMarconi et al.l |2005| ) . 



Three approaches are adopted to calibrate the Cepheid PL relations: cluster distances, 
trigonometric parallaxes and Baade-Wesselink (BW) method. Parallaxes based on Cepheids 
i n open clusters and a ssocia tions have be e n wid ely adopted to calibrate the PL relations 
( jTurner fc Burkd |2002| ) , but iFouque et al.l (120071 ) found a sys tematic difference wit h other 
methods. Accu r ate C epheid parallaxes have been provided by Benedict et al.l (120021 ) and by 
Benedict et al.l ( 120071 ) using the Fine Guide Sensor on board the Hubble Space Telescope 
(HST). They measured the distance of several Galactic Cepheids with average relative errors 
of ^ 8%. The BW method appea r s to be a robust approach, but some assumptions need 
to be investigated (lGautschylll987l : iBono et al.l 1 19941 : iGieren et al.ll2005l ). The value of the 
projection factor (p), i.e. the pararneter adopted to transform radial into pulsation velocity, 
is still controversial (ISabbev et al.l Il995l: iMerand et al.l l2006l ) . The typical value used in 
the literature is p = 1.36 (JBurki et al.lll982l ). but period-dependent values are also adopted 



(iGieren et al. 



by 



1993 



Merand et a 



( INardetto et al 



2004h 



20051). However, the p-factor of S Cephei, p = 1.27 ± 0.06, measured 



(120051) . using the HST parallax, agrees quite well with the predicted value 



Although, sev eral Cepheids are b inaries dynamical mass estimates are only available for 



a handful objects (JEvans et a 



20051). Current mass estimates rely ei t her on evolutionary 



or on pulsation prescriptions (JD'Cruz et al.ll2000l : iKeller &: Woodll2002l : ICaputo et al.l 1200 



and we are still facing the long-standing problem of the "Cepheid mass discrep ancv" (ICox 



2001 



1980 ) . Evolutionary mas s es are systeniatical ly larger than pulsation masses (JBono et al. 



Beauheu et al.ll200ll : ICaputo et al.ll2005r ). 



Theoretical fits of Ce pheid light curves provide robust constraints on their in trinsic 
parameters and distances (IWood et al.lll997l : iBono et al.ll2002l : Keller fc Woodll2006l ). This 
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approach, up to now, has only been applied to Bump Cepheids. We present, for the first 
time, a simultaneous fit of optical/NIR light, radial velocity, and angular diameter variations 
of the prototype S Cephei. We selected this object because accurate measurements of light, 
radius, velocity curves, and of t he p-factor are availab le and its geometric distance is known 
with an accuracy of 4% f273+j^ IJenedict et al.ll2002[ ). 



2. Empirical and theoretical framework 



The optic al (V, I) and the near -infrared (NIR) K-hand light curve were collected by iKiss 



(Il998l ) and by iBarnes et al.l (Il997l ). We selected these light curves because both the shape 
and the amplitude in the optical bands depend on surface temperature and radius variation, 
while in i^-band the d ependence on temper ature is significantly reduced. The radial velocity 
data wer e collected by Bersier et al.l ( 19941 1^1. while the angular diameter and the projection 
factor by iMerand et al.l ( 20051 ). The observables for 5 Cephei are summarized in Table 1. 



The adopted theoretical framework is based on a nonlinear radial puls ation code, includ- 



ing the nonlocal and time-dependent treatrn ent of turbulent convection ( (jStellingwerilll982 
Bono &: Stellingwerj 1 19941 : iBono et al.l Il999l ) . The system of nonlinear equations is closed 
using a free parameter, ami, that is proportional to the mixing length parameter. Changes 
in the ami parameter affect, as expected, both the limit cycle stabili t y (pulsation ampli- 



tudes ) and the topology of the instability strip (JFiorentino et al.l 120071 : iDi Criscienzo et al. 



20041). Similar app r oaches for the treatment of co nvective transpo r t hav e been developed by 



Feuchtinged (119991 ): iBuchler fc Szabol (120071 ). and lOliver fc Wood! (120051 ). 



In order to provide a detailed fit of 6 Cephei pulsation properties we adopted three 
different metal abundanc es, ranging from Z = 0.01 to Z = 0.02, to account for current spec- 
troscopic measurements (JFry &: Carneylll997l ) that indicate a metal abundance close to solar 
{[Fe/H] = —0.01 ± 0.06). Recent spectroscopic measurements of solar heavy element abun- 
dances show a significant decrease when compared with previous estima tes, and indeed the 
ratio between h eavy-eleraents and hydrogen , Z/X, decreased from 0.(3 23 (JGrevesse fc Sauval 
1998h to 0.0177 JLoddersl[2003l ) or to 0.0171 JAsplund et al.ll2004l . [2005h . This implies that the 
surface heavy element abundance of the Sun decreased from Z = 0.017 to Z = 0.0126/0.0133. 
The surface He abundance of the Sun according to helioseismological estimates, combine d 
with these latest solar metallicity estimates, is F = 0.2485 ± 0.0034 (JBasu fc Antiall2004j ). 
Howev er, the new solar photospheric abundan ces are at odds with helioseismology measure- 
ments (JBahcall et al.ll2005l : JGuzik et al.ll2005l ) and, at present, solar models are constructed 



^See the database |http://crocus.physics.mcmaster.ca/Cepheid/ 
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by adopting protosolar helium abundances ranging from Yq ~ 0.26 to Yq ~ 0.28 and heavy 
metal abundances ranging from Zq = 0.015 to Zq = 0.02. According to these evidence we 
adopted Z = 0.02 and Y = 0.28 (see Table 2). Once the chemical composition was fixed, we 
constructed a sequence of isoperiodic pulsation models, by adopting the pulsation relation 
(the relatio n connecting the p eriod to stellar mass, luminosity, and effective temperature) 
provided by lBono et al.l (120001 ). We selected the fundamental models with periods within 2% 



of the observed period (P=5.3663 days, Fernie 199q3) and pulsationally unstable ac c ordiri g 
to the analytical relations for the instability strip boundaries provided by lBono et al.l (120001 ). 
For each mass value, we fixed the l ower limit to the in put luminosity according to the Mass- 
Luminosity (ML) relation given by JBono et al.l (120001 ) and the upper limit according to mild 
overshooting p redictions, i.e. by adding to the quoted value a fixed luminosity excess i.e. 
AlogL = 0.25 (IChiosi et al.lll993l ). The bolometric light curves were transfo rmed into the 
UBVI JK photometric bands by adopting the model atmospheres provided by ICastelli et al. 
(1l997h . 



3. Results and discussion 

We constructed a sequence of models at fixed chemical composition (Z = 0.02, Y = 
0.28) and by assuming steps in mass and in effective temperature of AM = O.5M0 and 
ATe = lOOK. The comparison between theory and observations indicate that a model 
with M = 5.5Mq, and Te = 5700K accounts for the shape of the light curves and for the 
secondary bump located on the decreasing branch of the radial velocity curve (0 ~ 0.2). To 
further improve the agreement between predicted and empirical observables, we constructed 
a new sequence of models by adopting smaller steps in mass and in effective temperature, 
namely AM = 0.1MQ,ATe = 25i^. For every model in this sequence, the mixing length 
parameter was changed from ami = 0.3 to ami = 0.4. Note that an increase in ami causes 
a decrease in the pulsation driving, and in turn a decrease in the limit cycle amplitudes. 
Data plotted in the left panels of Fig. 1 show that the best fit model {ami = 0.33) accounts 
for the optical and NIR light curves, but both the radial velocity {a{AVr) = 4.2) and the 
radius (o"(Ai?) = 3.3) amplitudes are systematically smaller than the observed ones. This 
discrepancy can only be removed by decreasing the p-factor from 1.27 to the barely plausible 
value of p ~ 1. To further constrain the impact of the adopted chemical composition on 
the limit cycle stability, we constructed a new sequence of models by using the same metal 
abundance and a lower He content {Y = 0.26). The new best fit model (see Table 2 and the 
middle panels of Fig. 1) accounts for optical and NIR light curves and better reproduces 



^http://www. astro. utoronto.ca/DDO/research/cepheids/ 
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radial velocity (cr(AV^) = 1.6) and radius {a{AR) = 2.3) variations. The predicted \^-band 
light curve presents a mild discrepancy along the rising branch. However, the discrepancy on 
luminosity amplitudes is smaller than OAa. In order to verify the occurrence of a degeneracy 
in the input parameters, we also computed a new series of models by using a lower He content, 
(Y = 0.24, see Table 2 and right panels of Fig. 1). Interestingly enough, velocity and radius 
amplitudes predicted by new best fit-model are at least 2.5a smaller than observed. Once 
again, a reasonable fit between theory and observations would require a barely plausible 
p-factor {p = 1.03). 

In order to test how the fits depend on metallicity, we computed a new series of models 
by adopting Z = 0.015, and He content oi Y = 0.27 (see Table 2). Data plotted in the left 
panels of Fig. 2 show that the new best fit model agrees reasonably well with the observations. 
However, the discrepancy in the radial velocity is 1.6a and in the radius amplitude is 2.5a. 
This finding indicates that a decrease in the metal content improves the agreement between 
theory and observations. Therefore, we constructed two new sequences of pulsation models 
by adopting a more metal-poor chemical composition {Z = 0.01) and two different He 
contents {Y = 0.26,0.24). The new best fit models (see Table 2, middle and right panels of 
Fig. 2) account quite well for luminosity, radial velocity, and radius variations. A marginal 
difference between theory and observations is still present. However, the discrepancy is on 
average smaller than la for Y = 0.26 and systematically larger for Y = 0.24 (see Table 2). 
The comparison between theory and observations brings forward several findings, i) the more 
metal-poor {Z = 0.01) best fit models account for the observed light and velocity amplitudes 
at a level better than 0.8cr and the average discrepancy over the entire cycle is < 0.2S (see 
Table 2). ii) The fit of the radial velocity curve is very sensitive to the input parameters. 
The small "bump" located along the decreasing branch (0 = 0.2) mainly depends on the 
mass value, and it is independent of the chemical composition. Moreover, the radial velocity 
amplitude strongly depends on the He content. For example, an increase/decrease of ~ 8% in 
He content causes, at Z=0.02, an increase on the amplitude discrepancy of at least 1.4cr. The 
more metal-rich best fit models (Y=0.28, 0.25) to properly fit the velocity amplitude would 
require jo-factor valu es that are at least ~ 20% smaller than the empirical measurement by 



Merand et al.l (120061 ). Therefore, they do not pass the empirical validation, in) Predicted 
radius variations are at least 1.7a smaller than observed while the average discrepancy over 
the entire curve is never smaller than 0.2S for all the adopted chemical compositions. The 
mismatch is mainly caused by the few points located across the phases of minimum radius. 
We do not have a firm explanation for this discrepancy. These phases are the most rapid 
ones along the pulsation cycle and the radial displacemen ts of the outermost l ayers might be 



affected either by shock waves or by velocity gradients (JMathias et al.ll2006l ). This critical 



point deserves further investigations, since radius measurements depend on the adopted 



-6- 



distance and on limb darkening correction, but they are affected neither by the interstellar 
extinction nor by the p-factor. Therefore, they can play a key role in constraining theoretical 
predictions, iv) The model with Z = 0.01 and Y = 0.26 provides the optimum fit to all 
observables (except for the only three magnitudes on the rising branch of the V light curve). 
This finding taken at face value indicates that the He content might be 8% smaller than 
typically assumed to construct Galactic Cepheid models. The metal content of 6 Cephei 
seems to be 0.1—0.2 dex more metal-poor than suggested by spectroscopy, but this constraint 
is less robust. Note that the p-factor adopted for this model {p = 1.28) agrees quite well 
with the measurement by Merand et al. (2006). v) The best fit models present a luminosity 
excess (see column 7 in Table 2) that ranges from ALoqL/Lq = 0.07 (Z=0.01, Y=0.26) to 
0.17 (Z=0.02 , Y=0.24). These values are between the luminosities predicted by a canonical 



ML relation (JBono et al.ll2000l ) and those predicted by mild core overshooting ML relation 



{ALogL = 0.25, Chiosi et al. 1993). 

As a further validation of the current approach, we also estimated the absolute distance 
to S Cephei using the best fit models that passed the empirical validations. In order to provide 
a plausible estimate of the errors affecting distance determinations we assumed the largest 
uncertainties on the adopted input parameters: AM = O.IMq, AT = 50K, ALogL = 0.02, 
AR = 0.5Rq. By accounting for the entire error budget the intrinsic uncertainties on current 
distance determinations is about 3% for photometric distances and ~ 5 pc for the distances 
derived from the fit of radius variations. Distance estimates based on the apparent distance 
moduli prov ided by the fit of bo th optical and NIR light curves and on the reddening law 



provided by lCardelli et al.l (119891 ) are listed in columns 6 and 7 of Table 3. The comparison 
with the HST parallax shows that distances based on V, /-bands agree at the 0.2o" level, 
while those based on the V, /^-bands are within 0.8cr level. Errors either in the photometric 
zero-points or in the adopted color-temperature transformations or in the extinction law 
account for the difference in the distance determinations based on light curves. Mo reover, 



the occurrence of a circumstellar envelope (CSE) around 6 Cephei (JMerand et al.ll2006l ) could 
also affect the reddening corrections due to the different relative distribution of the absorbing 
grains around the star, compared to typical interstellar medium distribution. Independent 
dis tance determinations follow from comparing the observed angular diameter variations 



by iMerand et al.l (120051 ) with the radius curves from four different models. The distances 
(Table 3, column 8) agree well with each other and are on average only 5% (1.3o") larger 
than the trigonometric distance. The comparison of observed and model radii is in principle 
a powerful method to derive Cepheid distances, since it is independent of the assumption on 
the p-factor and on interstellar extinction. Note that an increase in the distance would also 
imply an increase in the p-factor, since at fixed radial velocity amplitude, p oc A9 • d, where 
AG is the angular diameter amplitude. 
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Finally we mention that the mean radii predicted by the best fit models match within 
lo" the empirical estimate based on the angular diameter provided by ME06 {R = 43±2Rq), 
by assuming the HST parallax. 

It is a pleasure to thank A. Merand for several interesting discussions on interferometric 
measurements. We also acknowledge an anonymous referee for a detailed and helpful report. 
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Table 1. Direct and derived observables of 6 Cephei. 


mv'^ 


mi'' rriK'' Ay'' Ai'' Ak'' AVr"" R"" AR^ ^J.^ p^ 


4.02 


3.23 2.34 0.84 0.53 0.20 47.8 43 ± 2 5.4 ±0.2 274 ± 11 1.27 ±0.06 



'^A pparent V, /, iC-band mean magnitudes a nd luminosity ampli tudes (|Kisslll998l : iBarnes et al.l 
1997^■ radial velo city amplitude ( km/sec. iBersier et al.l Il994h . mean radius (solar units , 
Merand et al.|2 006|). radius a mplitude iJMerand et al.l2006l) . true distance (pc. [Benedict et al.l2002h . 
p-factor JMerand et ai]|200dl . 



Table 2. Adopted intrinsic parameters of the best fit models and discrepancies between 

theory and observations. 



Y'a 



M" 



AL'' 



av^ 



err 



TK^ 



avr 



OR- 



0.02 



0.28 5.5 3.31 5800 45.06 0.11 0.33 O.i 



0.04 0.03 0.01 



0.4 



0.3 



0.9 0.3 0.3 4.2 



3.3 



0.02 


0.26 


5.5 


3.32 


5800 


45.07 


0.15 


0.40 


1.15 


0.2 


0.07 


0.008 


0.1 


0.2 


0.2 


0.4 


0.3 


1.6 


2.3 


0.02 


0.24 


5.5 


3.29 


5700 


45.32 


0.17 


0.38 


1.03 


0.08 


0.05 


0.01 


0.2 


0.3 


0.1 


0.5 


0.0 


3.0 


2.8 


0.015 


0.27 


5.5 


3.31 


5750 


45.29 


0.08 


0.40 


1.18 


0.2 


0.06 


0.009 


0.1 


0.2 


0.5 


0.7 


0.3 


1.6 


2.1 


0.01 


0.26 


5.5 


3.33 


5800 


46.07 


0.07 


0.43 


1.28 


0.2 


0.07 


0.01 


0.07 


0.2 


0.6 


0.8 


0.5 


0.4 


1.7 


0.01 


0.24 


5.5 


3.33 


5800 


46.08 


0.11 


0.44 


1.21 


0.1 


0.05 


0.009 


0.1 


0.2 


0.5 


0.7 


0.5 


1.0 


2.0 



''Metal (Z) and He (Y) abundance by mass, stellar mass {M/Mq), logarithmic l umin osity (logL/L0), effective temperature (K), radius 
(R/Rq), luminosity excess (dex) respect to canonical ML relation from lBono et al.l 1120001 ). convection efficiency parameter (omi), projection 
factor required to fit the observed amplitudes. 

''Standard deviations of the fit between the predicted curves and observations over the entire cycle. 

'^Standard deviations between observed and predicted amplitudes. 



Table 3. Distance and absolute magnitudine determinations based on different methods. 



Z'^ 


ya 


MV^' 


W" 


Mif'' 


di" 


d2^ 


da" 


Mv^ 


Mi^ 


A/if" 


0.02 


0.26 


7.50 


7.37 


7.28 


271 


282 


286 


-3.51 


-4.19 


-4.97 


0.015 


0.27 


7.46 


7.35 


7.29 


274 


284 


288 


-3.47 


-4.16 


-4.97 


0.01 


0.26 


7.52 


7.40 


7.32 


277 


288 


291 


-3.55 


-4.23 


-5.01 


0.01 


0.24 


7.52 


7.40 


7.32 


277 


288 


292 


-3.54 


-4.21 


-5.01 



"Metal (Z) and He (Y) abundance by mass; V, I, K apparent distance moduli 
(mag); true distances (pc) based on the Cardelli et al. (1989) extinction law; 
E{B-V) = (/iy-A'/)/l-22 (di), E{B-V) = (/^^ -Mk)/2.74 (^2); true distances 
(pc) based on angular diameter variations (ds); V,I,K absolute magnitudines.) 
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' ' ' I ' ' ' I ' ' ' I ' ' ' I ' ' ' 
Z = 0.02 Y=0.28 
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Fig. 1. — Comparison between best fit models at different chemical compositions and obser- 
vations. From top to bottom V, I, K light curves, radial velocity (km/sec), and radius curve 
(solar units). 
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' ' ' I ' ' ' I ' ' ' I ' ' ' I ' ' ' 

Z = 0.015 Y=0.27 
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Fig. 2. — Same as Fig. 1, but for more metal-poor best fit models. 



